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Abstract. We calculated the ellipticity of the deformed 
star due to the rotation or magnetic field. These two ef- 
fects are compared to each other within general relativity. 
It turned out that the magnetic distortion is important 
for recently observed candidates of magnetars, while the 
magnetic effect can be neglected for well-known typical 
pulsars. 
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1. Introduction 

Observations of pulsars have been accumulated since the 
first discovery by Hewish et 



al. (1968), and several new 
types of pulsars appeared with great surprise. These ob- 
servations have partially revealed the structure and evo- 
lution of rotating neutron stars. Their rotation periods 
range from 1.5 ms to several seconds. The surface mag- 
netic fields range from 10 8 to 10 15 G. The upper limit was 
recently raised by a factor 10 3 by the discovery of the 
soft gamma-ray repeaters (SGRs) and anomalous X-ray 
pulsars (AXPs) (se e e.g. Kouveliotou et al. 1998 , 1999[ 
Mereghetti & Stella 1995 ). The new class of pulsars with 
the very strong magnetic field (B ~ 10 14 -10 15 G) are most 
likely candidates of magnetars (e.g. Thompson & Dun- 
can 1996) and may be worth studying further. In the fu- 



ture, we may find a more extreme case, that is, a rapidly 
rotating relativistic star with a strong magnetic field. 

Most stars have spherically symmetric structure. They 
are however deformed due to the rapid rotation and the 
strong magnetic field. It is well known that both effects 
produce a flattening equilibrium star. We will examine 
them within general relativity. These effects are assumed 
to be small and treated as perturbations to spherically 
symmetric stars. The rotational axis of pulsars does not 
in general coincide with the axis of the dipole magnetic 



field. The relativistic treatments for the case is a com- 
plicated task, because the situation is not stationary. In 
this paper, however, we assume that the rotational effect 
decouples from the magnetic effect. Thus we consider the 
deformation arising from each perturbation separately and 
estimate the ellipticity. The estimate is important to judge 
which effect dominates in the rotating magnetized stars, 
whose rotation rate and the magnetic field are in a wide 
range. Our treatment is beyond the classical estimate in 
the Newtonian gravity, and give better comparison of the 
rotational effect and the magnetic effect on star deforma- 
tion for various pulsars. When either of them is huge, our 
estimate breaks down, and sophisticated numerical codes 
are required (see e .g. Bocquet et al. 1995 Bonazzola & 
Gourgoulhon 1996 ). 

The paper is organized as follows. In Sect. 2, we briefly 
review deformation of stars due to the rotation (Chan- 



drasekhar [L933|; Chandrasekhar & Roberts 

Fermi 

al 



magnetic field (Chandrasekhar 
Gal'tsov et 



1954 



1984 



1963 ) and the 
1953; Ferraro 
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Gal'tsov & Tsvetkov jl984| ) 
within Newtonian gravity. The quadrupole deformation 
can be evaluated by the ellipticity of the equilibrium 
shape. In Sect. 3, we also calculate the ellipticity based 
on the general relativistic perturbation theory (see also 
Hartle ( |1967D and Chandrasekhar fc Mi ller ( ]1974| ) for the 
rotational cases and Konno et al. ( |1999| ) for the magnetic 
cases). The ellipticity can be summarized in the same form 
as the Newtonian cases, but with different numerical fac- 
tors. In Sect. 4, using the ellipticity, we compare the rota- 
tional effect with the magnetic effect on star deformation 
numerically. In this comparison, we keep the parameter 
range of known pulsars in mind. Finally, we give conclud- 
ing remarks in Sect. 5. Throughout the paper, we use the 
units in which c = G = 1. 

2. Simple estimate of deformation 

Quadrupole deformation of an equilibrium body is char- 
acterized by the ellipticity e, which is defined by 
equatorial radius — polar radius 



mean radius 



(1) 



2 
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For the gravitational equilibrium with uniform rotation, 
the value is essentially related to the ratio of the rotational 
energy to the gravitational energy. In the slow rotation of 
a homogeneous star, we have a well-known result (see e.g. 



Chandrasekhar 1969): 



en 



4 M 



(2) 



where R, M and fl denote the radius, mass and angular 
velocity, respectively. For other stellar model, the numer- 
ical factor 5/4 should be replaced by an appropriate one. 
For example, the factor is 0.76 for the model with the poly- 
tropic equation of state (EOS) p — np^ n+1 ^ n with index 



n = 1 (see Table 1 in Chandrasekhar & Roberts 1963) 
We therefore generalize the expression with a dimension- 
less factor / as 



3r>2 



mi 



(3) 



and will discuss how the factor / depends on the stellar 
models. 

In a similar way, the effect of the magnetic stress is 
also expressed by the energy ratio of the magnetic field to 
the gravitational field. Introducing a dimensionless factor 
g, we have the ellipticity sb arising from magnetic field as 



£b = g 



APR 2 1 



(4) 



where /i is the magnetic dipole moment. The dimensionless 
factor g, in general, depends on both the magnetic field 
configurations and the EOS. For example, in the case of 
an incompressible fl uid b ody with a dipole magnetic field 
treated by Ferraro ( 1954 ), we derive g = 25/2. 



3. Relativistic calculation of deformation 

In this section, we will review quadrupole deformation due 
to the slow rotation and weak dipole magnetic field. The 
deformation can be expressed by the second-order quanti- 
ties with respect to the rotation rate or the magnetic field 
strength. In order to calculate the shape, we also have 
to calculate the space-time metric, which is axisymmetric 
stationary or static one. The line element can be written 
in the form 



ds 2 = - e v ^\l + 2 {h (r) + h 2 (r)P 2 (cos 9)}} dt 2 
2e A ( r 



,A(r) 



1 



{m (r) + m 2 (r)P 2 (cos 9)} 



+ r 2 (l + 2k 2 (r)P 2 {cos9)) 

x [d9 2 + r 2 sin 2 9(d<p - cu{r)dt) 2 ] , 



dr l 



(5) 



where P 2 is the Legendre's polynomial of degree 2, and 
(ho, h 2 , mo, m 2 , k 2 ) are the second-order quantities for the 
rotation rate or the magnetic field strength. The quantity 



to is the angular velocity acquired by an observer falling 
freely from infinity to a point r, and is equal to zero for 
the static magnetic field deformation. 

The stress-energy tensor of the perfect fluid body is 
described by 



(6) 



When we consider the magnetic field deformation, we fur- 
ther take into account the stress-energy tensor arising 
from the magnetic field, i.e., 



rp n 

J (cm) v 



(7) 



Solving the Einstein-Maxwell equations, we can obtain the 
second-order metric functions mentioned above. 

The ellipticity of the relativistic star can be calculated 
from the definition (Q) as 



3 j_ i. 

£ = -2 (T + k2 



(8) 



where £2 represents the displacement of quadrupole defor- 
mation. 

Since the displacement of the surface can be deter- 
mined by the hydrostatic equilibrium condition, the ellip- 
ticity of the slowly rotating star is expressed as (Chan- 
drasekhar & Miller |1974[ ) 



3 r 
£0 = —7/12 + ~;e- u <Sl 



h 2 . 

2 



(9) 
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Fig. 1. The dimensionless factor /, that is, the ellipticity of 
the rotational cases normalized by R 3 Q 2 /M, which is plotted 
against M/R. 

In order to compare it with the Newtonian results, we nor- 
malize the ellipticity in the same form as Eq. (0). There 
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are several possibilities of the normalization factors for 
M and R in the relativistic calculation. We use natural 
choices, i.e., gravitational mass for M and circumferential 
radius for R in this paper. This normalization is useful to 
extrapolate from the Newtonian results. Chandrasekhar & 



Miller (1974) used a different normalization, which causes 
a prominent peak (see Fig. 5 in their paper). We note that 
the peak is due merely to a less convenient choice of the 
normalization. The resultant dimensionless factor / is cal- 
culated for the polytropic EOS with n = 0, 0.5, 1, 1.5. 
Fig. |l| displays the variation of the dimensionless quantity 
/ with respect to the relativistic factor M/R. This figure 
shows that the correct relativistic calculations give smaller 
values of the ellipticity than those of the Newtonian cal- 
culations with fixed R 3 tt 2 /M. The factor in the typical 
relativistic case decreases down to 0.7 of the Newtonian 
case. 



variation of the ellipticity with respect to the relativis- 
tic factor M/R. In these calculations, we have used the 
current distribution 



J* oc po +po, 



(11) 



which is a simplest case derived by the integrability con- 
dition (see e.g. Ferraro 1954 ) an d give s the direct general- 
relativistic extension of Ferraro ( [1954 ) . Using the normal- 
ization, the residual factor g is almost independent of the 
relativistic factor M/R. 

4. The comparison 

In order to compare the rotational effect with the magnetic 
effect on star deformation, we now consider the ratio of en 
to £ B , 
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Fig. 2. The dimensionless factor g, that is, the ellipticity of 
the magnetic case normalized by p 2 / (M 2 R 2 ) , which is plotted 
against M/R. (Note that the normalization of this figure is 
different from that of Konno et al. (1999).) 



As for the weakly magnetized star, the Lorentz force 
plays a role on the equilibrium. The rotational term Q — lu 
is replaced by the magnetic term in the ellipticity. From 
the hydrostatic condition, we have (Konno et al. 1999) 



3 , 2r 
SB = — jfi2 + —, : r 



J* A* 



(10) 



where the subscript '0' denotes the background quantities, 
and and A^ are the ^-components of the 4-current and 
4-potential respectively. The ellipticity can also be writ- 
ten in the form of Eq. (Q). We also use the gravitational 
mass M and circumferential radius R. Fig. || displays the 



f_ 

.9 



R 3 n 2 

M 



\M 2 R< 



(12) 



0.10 



0.08 



§ 0.06 



0.04 



0.02 



n=0 ""-^ 




n=0.5 




n=Td — 




n=l.5 





0.00 



0.10 
M/R 



0.20 



Fig. 3. The ratio of the two dimensionless factors, f/g, plotted 
against M/R. 



First, we investigate the dependence of the ratio of the 
two dimensionless factors, f/g, on the relativistic factor 
M/R. Fig. H displays this ratio. From this figure, we can 
see that the true relativistic calculations with M/R ~ 0.2 
give smaller values of the ratio than those obtained by 
the Newtonian calculations. This fact means that the ap- 
proach using the Newtonian gravity overestimates the ro- 
tational effect. The curves f/g are approximately repro- 
duced within 10% if we use 



/ _ 1 - 1AM /R 
~g 10 + 8n 



(13) 
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Next, we consider the comparison including other pa- 
rameters f2 and fi of stars. We use the stellar model with 
M = 1.4Mq and B = 10km. Fig. || displays a critical line 
on which £q = £b and the two regions divided by this line 
in B-fl space, where B denotes the typical magnetic held 
strength on the surface, which is defined by B — /i/R 3 . 
We have plotted only one representative line of n — 1. We 
can also derive very close results for other indices. The 
critical line, in general, can be written from Eq. ( |l2| ) as 



B[G] » 5.3 x 10 1 



If ( M/lAMg 
g V i?/10km 



1/2 



Vl [sec 



(14) 



where it is useful to use the fitting formula Eq. (|13|) for 
f/g. In the region I, the magnetic effect dominates the 
rotational effect, i.e., £b > En, whereas in the region II 
vice versa, i.e., £r> > Eb- From this figure, we find that 
objects having mag netic held strength B ~ 10 14 -10 15 G 
and period T ~ 1 sec such as SGRs and AXPs belong 
to the region I. Thus, the magnetic effect overwhelms the 
rotational effect for such observed candidates of magne- 
tars. On the other hand, well-known typical pulsars with 
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Fig. 4. The critical line on which £n = ea in B-£l space for 
M/R = 0.2. This line divides the B-Q space into two regions. 
In the region I, the magnetic effect dominates the rotational 
effect, while in the region II vice versa. 



magnetic held strength B ~ 10 11 -10 13 G a nd the period 
T ~ 10 -1 -1 sec (see e.g. Taylor et al. 1993 ) obviously be- 
long to the region II. Millisecond pulsars also belong to 
the region II. The magnetic deformation is neglected. 



5. Concluding remarks 

The new classes of objects, which are candidates of mag- 
netars, have inspired us to investigate the relation between 



the rotational effect and the magnetic effect on deforma- 
tion of stars. We have briefly reviewed the quadrupole 
deformation due to the rotation and that due to the mag- 
netic field based on previous studies, and compared the ro- 
tational effect with the magnetic effect for various pulsars 
reported observationally. From our investigation, we have 
found that the new classes of objects such as SGRs and 
AXPs belong to the region in which the magnetic effect 
dominates the rotational effect, while well-known typical 
pulsars with magnetic field strength 10 11 -10 13 G and mil- 
lisecond pulsars vice versa. Thus, the critical line on which 
the ellipticity arising from the rotation equals to that aris- 
ing from the magnetic field divides the new classes and the 
well-known pulsars. Once we obtain the parameters O, B, 
R and M of a pulsar, we can see whether the magnetic 
effect is dominant or the rotational effect is dominant us- 
ing Eqs. (|l3|) and ( |l4| ) and Fig. ||[ The deformation due to 
the magnetic field may come into play in the spin-down 
of the magnetars. The spin-down history of the AXPs is 
rather irregular. Recently, Melatos (1999) ascribed the ir- 
regularity to the radiative precession produced by the de- 
formation. At present, the fit to the observational data 
is not so good, but will be improved by detailed models. 
Theoretical models for the deformed stars will be required 
there. 

As discussed by Bonazzola & Gourgoulhon (1996), 
the non-axisymmetric distortion is also important for the 
gravitational emission. The deformation can be calculated 
from the magnetic field, which is estimated from the ob- 
served pulsar period and the period derivative. The in- 
ferred amplitudes of the gravitational waves are too small 
for the present known pulsars. In the future, we might 
hnd more extreme case such as the early stage of rapidly 
rotating magnetars. 
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